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On " The process of collective de-excitation of atoms in a metastable level into emission mode of a single 
^ photon plus a neutrino pair, called radiative emission of neutrino pair (RENP), is sensitive to the 
absolute neutrino mass scale, to the neutrino mass hierarchy and to the nature (Dirac or Majorana) 
of massive neutrinos. We investigate how the indicated neutrino mass and mixing observables can 
be determined from the measurement of the corresponding continuous photon spectrum taking the 
example of a transition between specific levels of the Yb atom. The possibility of determining the 
nature of massive neutrinos and, if neutrinos are Majorana fermions, of obtaining information about 
the Majorana phases in the neutrino mixing matrix, is analised in the cases of normal hierarchical, 
inverted hierarchical and quasi-degenerate types of neutrino mass spectrum. We find, in particular, 
how critically the sensitivity to the nature of massive neutrinos depends on the atomic level energy 
difference relevant in the RENP. 



X 



Abstract 



1 Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 Sofia, Bulgaria 



1 



1 Introduction 



Determining the absolute scale of neutrino masses, the type of neutrino mass pattern, which can be 
either the normal or the inverted ordering (NO or 10), the nature (Dirac or Majorana) of massive 
neutrinos, and getting information about the Dirac and Majorana CP violation phases in the neutrino 
mixing matrix, are the most pressing and challenging problems of the future research in the field of 
neutrino physics (see, e.g., [I]). At present we have compelling evidence for existence of mixing of 
three massive neutrinos z/j, i — 1,2, 3, in the weak charged lepton current (see, e.g., [2]). The masses 
vrti > of the three light neutrinos V{ do not exceed a value approximately 1 eV, rrii ^ 1 eV. The three 
neutrino mixing scheme is described (to a good approximation) by the Pontecorvo, Maki, Nakagawa, 
Sakata (PMNS) 3x3 unitary mixing matrix, C/pmns- m ^ ne widely used standard parametrisation 
[I], ^pmns is expressed in terms of the solar, atmospheric and reactor neutrino mixing angles 6* 12 , 
6*23 and 6*13, respectively, and one Dirac (5), and two Majorana [3111] (a and (3) CP violation (CPV) 
phases. In this parametrisation, the elements of the first row of the PMNS matrix, U e i, i = 1,2,3, 
which play important role in our further discussion, are given by 

Uei = ci 2 c 13 , U e2 = s 12 c 13 e ia , U e3 = s 13 e 1 ^ , (1) 

where we have used the standard notation Cy = cos^-, s^- = sin^- with < Qij < tt/2, < 6 < 2n 
and, in the case of interest for our analysis □, < a, (3 < ir, (see, however, [5]). If CP invariance 
holds, we have 5 = 0, ir, and [6] a, (3 = 0, ir/2, ir. 

The neutrino oscillation data, accumulated over many years, allowed to determine the parame- 
ters which drive the solar and atmospheric neutrino oscillations, Awig = Am^, 6*12 and |Am^| = 
[Am^l = |Am| 2 |, 6*23, with a high precision (see, e.g., [2]). Furthermore, there were spectacular 
developments in the last year in what concerns the angle 9i 3 (see, e.g., ppj). They culminated in a 
high precision determination of sin 2 26* 13 in the Daya Bay experiment using the reactor v e [7]: 

sin 2 20i3 = 0.089 ±0.010 ±0.005. (2) 

Similarly, the RENO, Double Chooz, and T2K experiments reported, respectively, 4.9a, 3.1a and 
3.2er evidences for a non-zero value of #13 [8], compatible with the Daya Bay result. 

A global analysis of the latest neutrino oscillation data presented at the Neutrino 2012 Inter- 
national Conference [2] was performed in [9]. We give below the best fit values of Am^, sin 2 #i 2 , 
lAm 2 ^^! and sin 2 ^13, obtained in [9], which will be relevant for our further discussion: 

Am 2 ! = 7.54 x 10" 5 eV 2 , |Am 2 1(32) | = 2.47 (2.46) x 10~ 3 eV 2 , (3) 
sin 2 9 12 = 0.307, sin 2 6 Vi = 0.0241 (0.0244) , (4) 

where the values (the values in brackets) correspond to NO (IO) neutrino mass spectrum. We will 
neglect the small differences between the NO and IO values of lAm 2 ^^! and sin 2 #13 and will use 
|Am 2 -^ 32 )| = 2.47 x 10~ 3 eV 2 , sin 2 9 13 = 0.024 in our numerical analysis. 

After the successful measurement of #13, the determination of the absolute neutrino mass scale, 
of the type of the neutrino mass spectrum, of the nature of massive neutrinos, as well as getting 
information about the status of CP violation in the lepton sector, remain the highest priority goals 
of the research in neutrino physics. Establishing whether CP is conserved or not in the lepton sector 

2 We use the convention adopted in [T]. 

3 Note that the two Majorana phases «2i and a^i defined in [1] are twice the phases a and /3: (X21 = 2a, CK31 = 2/3. 
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is of fundamental importance, in particular, for making progress in the understanding of the origin 
of the matter-antimatter asymmetry of the Universe (see, e.g., [TU | [TT | fT2"]). 

Some time ago one of the present authors proposed to use atoms or molecules for systematic 
experimental determination of the neutrino mass matrix [131 EE]- Atoms have a definite advantage 
over conventional target of nuclei: their available energies are much closer to neutrino masses. The 
process proposed is cooperative de-excitation of atoms in a metastable state. For the single atom 
the process is |e) — > \g) + 7 + (z^ + Uj), i,j = 1, 2, 3, where i/j's are neutrino mass eigenstates. If z/j 
are Dirac fermions, (z/j + z/j) should be understood for i = j as (z/j + z/j), and as either (z/j + z/j) or 
[vj + z/j) when i 7^ j, z/j being the antineutrino with mass m;. If z/j are Majorana particles, we have 
z/j = vi and (z/j + z/j) are the Majorana neutrinos with masses m; and rrij. 

The proposed experimental method is to measure, under irradiation of two counter-propagating 
trigger lasers, the continuous photon (7) energy spectrum below each of the six thresholds ojy cor- 
responding to the production of the six different pairs of neutrinos, V\V\, UiU 2 ,..., z/ 3 z/ 3 : u < Uij, u 
being the photon energy, and 



e eg (im + rrij) 2 . . 
Uij = Uji = — , i,j = 1,2,3, mi, 2i3 >0, (5) 

where e e5 is the energy difference between the two relevant atomic levels. 

The process occurs in the 3rd order (counting the four Fermi weak interaction as the 2nd order) 
of electroweak theory as a combined weak and QED process, as depicted in Fig. [TJ Its effective 
amplitude has the form of 

^G F V. • a^v\avi 
(g\d\p) ■ E ■ (p\S e \e) , (6) 

e pg - cu 

a ij = U*iU ej - ^ , (7) 

where U e i, i = 1,2, 3, are the elements of the first row of the neutrino mixing matrix Upmns, given in 
eq. (OQ). The atomic part of the probability amplitude involves three states |e), \g), \p), where the two 
states |e), \p), responsible for the neutrino pair emission, are connected by a magnetic dipole type 
operator, the electron spin S e . The \g) — \p) transition involves a stronger electric dipole operator 
d. From the point of selecting candidate atoms, El x Ml type transition must be chosen between the 
initial and the final states (|e) and \g))- The field E in eq. ([6]) is the one stored in the target by 
the counter-propagating fields. The formula has some similarity to the case of stimulated emission. 
By utilizing the accuracy of trigger laser one can decompose, in principle, all six photon energy 
thresholds at Uij, thereby resolving the neutrino mass eigenstates instead of the flavor eigenstates. 
The spectrum rise below each threshold u < depends, in particular, on \a,ij\ 2 and is sensitive 
to the type of the neutrino mass spectrum, to the nature of massive neutrinos, and, in the case of 
emission of two different Majorana neutrinos, to the Majorana CPV phases in the neutrino mixing 
matrix (see below on further details). 

The disadvantage of atomic targets is their smallness of rates which are very sensitive to available 
energy of order eV. This can be overcome by developing, with the aid of a trigger laser, macro- 
coherence of atomic polarization to which the relevant amplitude is proportional, as discussed in 
[T5| [T6] . The macroscopic polarization supported by trigger field gives rise to enhanced rate oc n 2 V, 
where n is the number density of excited atoms and V is the volume irradiated by the trigger laser. 
The proposed atomic process may be called radiative emission of neutrino pair, or RENP in short. 
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Figure 1: A— type atomic level for RENP |e) — > \g) + 7 + v{Uj with z/, a neutrino mass eigenstate. 
Dipole forbidden transition |e) — > \g) + 7 + 7 may also occur via weak Elx Ml couplings to \p). 



The estimated rate roughly of order mHz or a little less makes it feasible to plan realistic RENP 
experiments for a target number of order of the Avogadro number, within a small region of order 
1 ~ 10 2 cm 3 , if the rate enhancement works as expected. 

The new atomic process of RENP has a rich variety of neutrino phenomenology, since there 
are six independent thresholds for each target choice, having a strength proportional to different 
combinations of neutrino masses and mixing parameters. In the present work we shall correct the 
spectrum formula for the Majorana neutrino case given in [T4] and also extend the discussion of the 
atomic spin factor. 

In the numerical results presented here we show the sensitivity of the RENP related photon 
spectral shape to various observables; the absolute neutrino mass scale, the type of neutrino mass 
spectrum, the nature of massive of neutrinos and the Majorana CPV phases in the case of massive 
Majorana neutrinos. All these observables can be determined in one experiment, each observable 
with a different degree of difficulty, once the RENP process is experimentally established. For atomic 
energy available in the RENP process of the order of a fraction of eV, the observables of interest can 
be ranked in the order of increasing difficulty of their determination as follows: 

(1) The absolute neutrino mass scale, which can be fixed by, e.g., measuring the smallest photon 
energy threshold min(oOij) near which the RENP rate is maximal: mm(ujij) corresponds to the pro- 
duction of a pair of the heaviest neutrinos (max(mj) ^ 50 meV). 

(2) The neutrino mass hierarchy, i.e., distinguishing between the normal hierarchical (NH), inverted 
hierarchical (IH) and quasi-degenerate (QD) spectra, or a spectrum with partial hierarchy (see, e.g., 



(3) The nature (Dirac or Majorana) of massive neutrinos. 

(4) The measurement on the Majorana CPV phases if the massive neutrinos are Majorana particles. 
The last item is particularly challenging. The importance of getting information about the Ma- 
jorana CPV violation phases in the proposed RENP experiment stems, in particular, from the pos- 
sibility that these phases play a fundamental role in the generation of the baryon asymmetry of 
the Universe [IT]. The only other experiments which, in principle, might provide information about 
the Majorana CPV phases are the neutrinoless double beta ((P(3)ou-) decay experiments (see, e.g., 



The paper is organized as follows. In Section 2 the basic RENP spectral rate formula is given 
along with comments on how the Majorana vs Dirac distinction arises. We specialize to rates under 



DP)- 
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no magnetic field so that the experimental setup is simplest. In Section 3 we discuss the physics 
potential of a RENP experiment for measuring the absolute neutrino mass scale and determining 
the type of neutrino mass spectrum (or hierarchy) and the nature (Dirac or Majorana) of massive 
neutrinos. This is done on the examples of a candidate transition of Yb J = metastable state 
and of a hypothetical atom of scaled down energy of the transition in which the photon and the two 
neutrinos are emitted. Section 4 contains Conclusions. 

2 Photon Energy Spectrum in RENP 

When the target becomes macro-coherent by irradiation of trigger laser, RENP process conserves 
both the momentum and the energy which are shared by a photon and two emitted neutrinos resulting 
in the threshold relation (j5]) pj)] . The atomic recoil can be neglected to a good approximation. Since 
neutrinos are practically impossible to measure, one sums over neutrino momenta and helicities, and 
derives the single photon spectrum as a function of photon energy u. We think of experiments that 
do not apply magnetic field and neglect effects of atomic spin orientation. The neutrino helicity 
(denoted by h r ,r = 1,2) summation in the squared neutrino current j k = dijv\<JkVj gives bilinear 
terms of neutrino momenta (see [T3] and the discussion after eq. (1171) ): 



hi,h,2 



. mirrij ( (Im(a i7 )) 2 \ \ 1 , k , _ _ , 



The case 5m = 1 applies to Majorana neutrinos, 5m = corresponds to Dirac neutrinos. The term 
oc mirrij(l — 2(Im(ai ? )) 2 /|aj ? | 2 is similar to, and has the same physical origin as, the term oc MjMj 
in the production cross section of two different Majorana neutralinos Xi an d Xj with masses Mj and 
Mj in the process of e~ + e + — > Xi + Xj [19] - The term oc MjMj of interest determines, in particular, 
the threshold behavior of the indicated cross section. 

The subsequent neutrino momentum integration (with Ei = \Jpf + m 2 being the neutrino energy) 

/ ^j0^ 5Z{% + pl + p 2)6{ee9 - u - Ei - E2)K ^ = hj dVvK ^ ' (9) 

can be written as a second rank tensor of photon momentum, kikjG^ + 5ijk 2 G^ from rotational 
covariance. Two coefficient functions G^> are readily evaluated by taking the trace Yli=j an d a 
product with kikj and using the energy-momentum conservation. But their explicit forms are not 
necessary in subsequent computation. 

We now consider sum over magnetic quantum numbers of ElxMl amplitude squared: 



R= f dV v ^^-J2\J2(9 M 9 \d- E\pM p ) ■ (pM p \S e -l\eM e )\ 2 . (10) 

J e M, 



The field E is assumed to be oriented along the trigger axis taken parallel to 3— axis. Since there is 
no correlation of neutrino pair emission to the trigger axis, one may use the isotropy of space and 
replace (S e ■ k) n i(S e ■ k)i n i by (S e ) n i ■ (S e )i n 'k 2 /3. Using the isotropy, we define the atomic spin factor 
C ep (X) of X atom by 



2J e + 1 



(pM p \S e \eM e ) ■ (eM e \S e \ P M' p ) = 5 MpM ,C ep {X) . (11 
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This means that only the trace part of eq. (|S]), AK^/3, is relevant for the neutrino phase space 
integration. 

The result is summarized by separating the interference term relevant to the case of Majorana 
neutrinos v^. 



r 72 » = rvMUt) , r = /< P9 e9 (2j p + i)c ep , (12) 



3n 2 VG 2 F jp g € eg n 

2e 3 
pg 

= — y~] \ a ij\ 2 feM _ TXlj B?f ) , (13) 

{e pg - u) *-r? 

2 0M*£l \ a -U*U h 



(14) 



A *iM = — 7 7TT { ( e e 9 (e e9 - 2w) - (m, + m 3 ) 2 ) (e e9 (e e9 - 2u) - (m, - rrij) 2 ) } 1/2 , (15) 



J *iM = Qe es (e efl - 2w) + ^w 2 - ^w 2 A 2 (u;) - i(m 2 + m 2 ; 



6 e lg( e eg ~ 2w) 



2 

/ 2 2\ f 



l^eg , o 2\2 



(16) 



The term oc 5m rriimj appears only for the Majorana case. We shall define and discuss the dynamical 
dimensionless factor r) u (t) further below. The limit of massless neutrinos gives the spectral form, 

J( " ;m - = 0)= 12(e pg -^ ' (17) 

where the prefactor of J^.. ja^ | 2 = 3/4 is calculated using the unitarity of the neutrino mixing matrix. 
On the other hand, near the threshold these functions has the behavior oc y/uJij — cu. 

We will explain next the origin of the interference term for Majorana neutrinos. The two- 
component Majorana neutrino field can be decomposed in terms of plane wave modes as 

^ M (x, t) = J2 (u(p)e- iEit+i P%(p) + u c (p)e" 5, *- i * 5f &J(p)) , (18) 

i,p 

where the annihilation bi(p) and creation b\{p) operators appears as a conjugate pair of the same 
type of operator b in the expansion (the index % gives the %— th neutrino of mass rrii, and the helicity 
summation is suppressed for simplicity). The concrete form of the 2-component conjugate wave 
function u c oc ia^u* is given in [13]. A similar expansion can be written in terms of four component 
field if one takes into account the chiral projection (1 — 7s)/2 in the interaction. The Dirac case is 
different involving different type of operators bi(p) and d\(p): 

^ D (x, t) = J2 (u(p)e- iEit+ip %(p) + v[p)e iEit -^ s d\{p) \ . (19) 



Neutrino pair emission amplitude of modes ipi,jp2 contains two terms in the case of Majorana 
particle: 

b\b] (^(plKG^) - a iiU *(^)u c (pi)) , (20) 



and its rate involves 
1 
2 



' 'ay«*(piK(£) -a^ipzXfa)] 2 = \\a l3 \ 2 (|^(1, 2)| 2 + |^(2, 1)| 2 ) - K(o 2 .) (^(1, 2)^(2, 1) 



(21) 
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where the relation = a*- is used and ij}(l,2) = u*(pi)u c (p 2 ). The result of the helicity sum 
^2 (V'(l) 2)^(2, 1)*) is in [T3], which then gives the interference term oc Bff in the formula (114jl . 



We see from eqs. (112jl and (1131) that the overall decay rate is determined by the energy independent 
To, while the spectral information is in the dimensionless function I(oS). The rate Tq given here is 
obtained by replacing the field amplitude E of Eqj6] squared by e eg n, which is the atomic energy 
density stored in the upper level |e). 

The dynamical factor (t) is defined by a space integral of a product of macroscopic polarization 
squared times field strength, both in dimensionless units, 



QmL/2 J ^ 1 (e,a m t) 2 + r 2 (e,a m t) 2 

a m L/2 



W) = 7 I d £ 3 |eU,a m t)| . (22) 



Here ri ± ir 2 is the medium polarization normalized to the target number density. 

The dimensionless field strength |e(£,r)| 2 = \E(£ = a m x,r = a m t)\ 2 / (e eg n) is to be calculated 
using the evolution equation for field plus medium polarization in [16], where £ = a m x (a m = 
e eg n ge n/2 with n ge the off-diagonal coefficient of AC Stark shifts [T3]) is the atomic site position in 
dimensionless unit along the trigger laser direction (— L/2 < x < L/2 with L the target length), 
and r = a m t is the dimensionless time. The characteristic unit of length and time are a" 1 ~ 
lcm(n/10 21 cm -3 ) -1 and 40ps(n/10 21 cm" 3 )" 1 for Yb discussed below. We expect that r] u (t) in the 
formula given above is roughly of order unity or lessQ We shall have more comments on this at the 
end of this section. 

Note that what we calculate here is not the differential spectrum at each frequency, instead it 
is the spectral rate of number of events per unit time at each photon energy. Experiments for the 
same target atom are repeated at different frequencies u)\ < Wn in the NO case (or u\ < W33 in 
the 10 case) since it is irradiated by two trigger lasers of different frequencies of U{ (constrained by 
uj\ + uji = e eg ) from counter-propagating directions. 

As a standard reference target we take Yb atom and the following de-excitation path, 

Yb; \e) = (QsQp) 3 P , \g) = (6s 2 ) % , \p) = (QsQp) 3 P 1 . (23) 

Relevant atomic parameters are as follows |20j : 

e eg = 2.14349 eV , e pg = 2.23072 eV , j pg = 1.1 MHz . (24) 

The notation based on LS coupling is used for Yb electronic configuration, but this approximation 
must be treated with care, since there might be a sizable mixing based on jj coupling scheme. 
The relevant atomic spin factor C ep (Yb) is estimated, using the spin Casimir operator within an 
irreducible representation of LS coupling. Namely, 

( 3 P \S e \ 3 P 1 ,M) ■ ( 3 P 1 ,M\S e \ 3 P ) = ^fPol^Px.M) • ( 3 P 1 ,M\S e \ 3 P ) = |, (25) 

M 

since S e ■ S e — 2 for the spin triplet. This gives C ep (Yb)= 2/3 for the intermediate path chosen. 

We also considered another path, taking the intermediate state of Yb, X P\ with e pg = 3.10806 eV , 
Ipg = 0.176 GHz. Using a theoretical estimate of A-coefficient 4.6 x 10 -2 Hz for l P\ — > 3 P\ transition 



4 There is a weak dependence of the dynamical factor r) u {f) on the photon energy ui, since the field e in Eq l221 
a solution of the evolution equation, is obtained for the initial-boundary condition of frequency u> dependent trigger 
laser irradiation. 
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given in NIST [2U] and taking the estimated Lande g-factor [21], 3/2 for the 3 Pi case, we calculate 
the mixed fraction of jj coupling scheme in LS forbidden amplitude squared |( 1 Pi|S' e | 3 Pi)| 2 , to give 
C ep ~ 1 x 10" 4 . 

Summarising, the overall rate factor T is given by 

r ° = o 3 2J p + 1 C e P ~ 0.37 mHz - — y— - , 26 

2e pg \\) zl cm d 10^ cur 

where the number is valid for the Yb first excited state of J = 0. If one chooses the other intermediate 
path, 1 P 1 , the rate T is estimated to be of order, 1 x 10~ 3 mHz, a value much smaller than that 
of the 3 Pi path. The denominator factor l/(e pg — w) 2 is slightly larger for the 3 Pi path, too. We 
consider the intermediate 3 Pi path alone in the following. 

The high degree of sensitivity to the target number density n seems to suggest that solid envi- 
ronment is the best choice. But de-coherence in solids is fast, usually sub-pico seconds, and one has 
to verify how efficient coherence development is achieved in the chosen target. 

Finally, we discuss a stationary value of time independent r] u (t) ( 122]) some time after trigger 
irradiation. The stationary value may arise when many soliton pairs of absorber-emitter [TH] are 
created, since the target in this stage is expected not to emit photons of PSR origin (due to the 
macro-coherent |e) — > \g) + 77), or emits very little only at target ends, picking up an exponentially 
small leakage tail. This is due to the stability of solitons against two photon emission. Thus the 
PSR background is essentially negligible. According to [22], the r)u(t) integral ( 122]) is time dependent 
in general. Its stationary standard reference value may be obtained by taking the field from a 
single created soliton. This quantity depends on target parameters such as a m and relaxation times. 
Moreover, a complication arises, since many solitons may be created within the target, and the 
number of created solitons should be multiplied in the rate. This is a dynamical question that has 
to be addressed separately. In the following sections we compute spectral rates, assuming rj u (t) = 1. 



3 Sensitivity of the Spectral Rate to Neutrino Mass Ob- 
servables and the Nature of Massive Neutrinos 

We will discuss in what follows the potential of an RENP experiment to get information about the 
absolute neutrino mass scale, the type of the neutrino mass spectrum and the nature of massive 
neutrinos. We begin by recalling that the existing data do not allow one to determine the sign of 
Am\ = Am^jj and in the case of 3-neutrino mixing, the two possible signs of Am^^ corresponding 
to two types of neutrino mass spectrum. In the standard convention [1] the two spectra read: 
i) spectrum with normal ordering (NO): mi < m 2 < m 3 , Am\ = Am 2 ^ > 0, Am^ > 0, 7712(3) = 
(m 2 + Am^^j)^; ii) spectrum with inverted ordering (IO): m 3 < mi < m 2 , Am\ = Am 2 2 < 0, 

Am\ Y > 0, m 2 = (m§ + Arn^) 5 , uii = (m 2 + Am| 3 — Am^ji Depending on the values of the 
smallest neutrino mass, min(raj) = m , the neutrino mass spectrum can also be normal hierarchical 
(NH), inverted hierarchical (IH) and quasi-degenerate (QD): 

NH : mi < m 2 < m 3 , m 2 = (Am 2 ^ ^ 0.009eV, m 3 = (Am 2 ^ = 0.05eV, (27) 
IH : m 3 < mi < m 2 , mi i2 = |Am 2 2 |5 ^ 0.05 eV, (28) 
QD : mi = m 2 = m 3 = m , m) > | Am 2 1(32) | , m £ 0.10 eV . (29) 

All three types of spectrum are compatible with the existing constraints on the absolute scale of 
neutrino masses m,-. 



8 



Table 1: The quantity |a.jj| 2 = \U* e JJ e j — 



\ a \2 _ |J r 2 _ 112 
|"11| — I c 12 c 13 2l 


|"12| — c 12 i 12 c 13 


1 n 1 1 ^ — ^2 
|«13| — c 12 6 13 c 13 


0.0311 


0.2027 


0.0162 


U 12 _ |„2 2 1 12 
l a 22 | — |Sl2 C 13 ~~ 2I 


^23 2 = Sl2 S 13 C 13 


\ n 12 _ 1 „2 1 12 
l fl 33| - |Si3 - 2I 


0.0405 


0.0072 


0.2266 



3.1 General features of the Spectral Rate 

The first thing to notice is that the rate of emission of a given pair of neutrinos (z/j + z/j) is suppressed, 
in particular, by the factor \a,ij\ 2 , independently of the nature of massive neutrinos. The expressions 
for the six different factors \dij\ 2 in terms of the sines and cosines of the mixing angles #12 and #13, as 
well as their values corresponding to the best fit values of sin 2 9 12 and sin 2 613 quoted in eq. d3J), are 
given in Table [TJ It follows from Tabled] that the least suppressed by the factor \aij\ 2 is the emission 
of the pairs (1/3 + z/3) and iv\ + v 2 ), while the most suppressed is the emission of (1/2 + z/3). The values 
of |dij| 2 given in Table [1] suggest that in order to be able to identify the emission of each of the six 
pairs of neutrinos, the photon spectrum (i.e., the RENP spectral rate) should be measured with a 
relative precision not worse than approximately 5 x 10~ 3 . 

As it follows from eqs. ffl3|) and (JHJ), the rate of emission of a pair of Majorana neutrinos with 
masses rrii and rrij differs from the rate of emission of a pair of Dirac neutrinos with the same masses 
by the interference term oc mirrijB^ . For % = j we have Bf- = 1, the interference term is negative 
and tends to suppress the neutrino emission rate. In the case of i ^ j, the factor Bf? , and thus the 
rate of emission of a pair of different Majorana neutrinos, depends on specific combinations of the 
Majorana and Dirac CPV phases of the neutrino mixing matrix: from eqs. (fl4"]l and ([!]) we get 

B\[ = cos 2a , = cos 2(/3 - 6) , B% = cos 2(a - f3 + 6) . (30) 

In contrast, the rate of emission of a pair of Dirac neutrinos does not depend on the CPV phases of the 
PMNS matrix. In the case of CP invariance we have a, /3 — 0, 7r/2, it, 5 = 0, it, and, correspondingly, 
Bf? = — 1 or +1, i ^ j. For Bf? = +1, the interference term tends to suppress the neutrino emission 
rate, while for BfJ = —1 it tends to increase it. If some of the three relevant (combinations of) 
CPV phases, say a, has a CP violating value, we would have — 1 < Bf£ < 1; if all three are CP 
violating, the inequality will be valid for each of the three factors Bfj: — 1 < Bfj < 1, i ^ j . Note, 
however, that the rates of emission of {v\ + u 3 ) and of (u 2 + ^3) are suppressed by |ai 3 | 2 = 0.016 and 
|a 2 3 1 2 = 0.007, respectively. Thus, studying the rate of emission of {v\ + ^2) seems the most favorable 
approach to get information about the Majorana phase a, provided the corresponding interference 
term oc rnxmiB^ is not suppressed by the smallness of the factor m\m2- The mass m\ can be very 
small or even zero in the case of NH neutrino mass spectrum, while for the IH spectrum we have 
m\m 2 ^ |Am3 2 | = 2.5 x 10 -3 eV 2 . We note that all three of the CPV phases in eq. (130]) enter into 
the expression for the (PP)ou— decay effective Majorana mass as their linear combination (see, e.g., 
HM): 

E 7-7-212 24, 244, 244 

rriiUJ = m 3 s 13 + m 2 s 12 c 13 + m lCl2 c 13 

i 

+2mim 2 s 2 2 c 2 2 c^ 3 cos(2a) + 2m 1 m 3 s 2 3 c 2 2 c 2 3 cos 2(/3 — 5) + 2m2m 3 s 2 3 s 2 2 c 2 3 cos 2(a — (3 + S) . (31) 
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Yb RENP global spectrum 




In the case of ra\ < m 2 < m 3 (NO spectrum), the ordering of the threshold energies at = Uji 
is the following: wy > cui 2 > ^22 > ^13 > ^23 > u 33 . For NH spectrum with negligible mi which can 
be set to zero, the factors (mi + mj) 2 = Kij in the expression (jSJ) for the threshold energy Uij are given 
by: k u = 0, «i2 = Am^, k 22 = 4Am| 1 , k 13 = Am^, k 2 3 = (-^/AmfJ" + y / Am^") 2 , k 33 = 4Amf 1 . It 
follows from eq.® and the expressions for that Wn, Ui 2 and u 22 are very close, u± 3 and W23 are 
somewhat more separated and the separation is the largest between u 22 and W13, and W23 and u 33 : 

- (cji 2 - w 22 ) = ^— Am 2 21 = 1.759 (8.794) x 10" 5 eV , (32) 



NH : 


0J n 


- Wl2 


NH : 


W13 


— <^23 


NH : 


^22 


- W13 


NH : 


<^23 


— W33 



1 



(2 y Amli y Amli + Am 2 21 ) = 0.219 (1.095) x 10" 3 eV , (33) 
(Am 2 31 - AAm 2 21 ) ^ 0.506 (2.529) x 10" 3 eV , (34) 



- — (3Am^ - 2 JAm 2 21 JAm^ - Am z 21 ) = 1.510 (7.548) x 10~ 3 eV , (35) 
2e e3 

where the numerical values correspond to Am 2l given in eq. ([3]) and e eg = 2.14349 ( numbers in 
parenthesis corresponding to the 1/5 of Yb value, namely 0.42870) eV. We get similar results in what 
concerns the separation between the different thresholds in the case of QD spectrum and Am 31 > 0: 

QD : uu - u 12 = u 12 - u 22 wig - w 2 3 = — Am 2 2l 3.518 (17.588) x 10" 5 eV , (36) 

QD : u 22 - u 13 ^ u 23 -u 33 Am 2 21 = — (Am 2 31 - 2Am 2 21 ) = 1.082 (5.410) x 10~ 3 eV(.37) 

(-eg ^eg 

For spectrum with inverted ordering, m 3 < mi < m 2 , the ordering of the threshold energies is 
different: u 33 > Ui 3 > u 23 > u u > u 12 > uj 22 . In the case of IH spectrum with negligible m 3 = 0, 
we have: k 33 = 0, K13 = Am 23 — Am 21 , k 23 = Am 23 , ku = 4(Am2 3 — Am 2l ), k± 2 = (\/ Am 23 + 
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Threshold behavior (Dirac case) 




1.064 1.066 1.068 1.070 

a) [eV] 



Figure 3: Photon energy spectrum from Yb 3 Pq — > 1 S'o in the threshold region for 3 different Dirac 
masses, mo = 2, 20, 50meV both in the NH(in solid) and IH(in dashed) cases. 



\J A77123 — Am^) 2 , k 2 2 = 4Am2 3 . Now not only uu, U12 and U22, but also W13 and U23, are very 
close, the corresponding differences being all ~ Am^/e^. The separation between the thresholds 
W33 and W13, and between U23 and uu, are considerably larger, being ~ Am^/fej. These results 
remain valid also in the case of QD spectrum and Am| 2 < 0. 

This precision is possible in our RENP experiments since the energy resolution in the spectrum 
is determined by accuracy of the trigger laser frequency, which is much better than 10 _5 eV. 

3.2 Neutrino Observables 

We will concentrate in what follows on the analysis of of the dimensionless spectral function I (u) 
which contains all the neutrino physics information of interest. 

In Fig. [2] we show the global features of the photon energy spectrum for the Yb 3 P 1 So 
transition in the case of massive Dirac neutrinos and NH and IH spectra. For m < 20 meV, all 
spectra (including those corresponding to massive Majorana neutrinos which are not plotted) look 
degenerate owing to the horizontal and vertical axes scales used to draw the figure. 

The Absolute Neutrino Mass. Much richer physics information is contained in the spectrum 
near the thresholds ojy. Figure [3] shows the Dirac neutrino spectra for three different sets of values 
of the neutrino masses (corresponding to the smallest mass m = 2, 20, 50 meV) and for both the 
NO (Ami^g^ > 0) and IO (Am^^j < 0) neutrino mass spectra. One sees that the locations of 
the thresholds corresponding to the three values of m (and that can be seen in the figure) differ 
substantially. This feature can be used to determine the absolute neutrino mass, including the 
smallest mass, as evident in differences of spectrum shapes for different masses of m , 2, 20, 50 meV 
in FigJ3J In particular, the smallest mass can be determined by locating the highest threshold (wu 
for NO and W33 for IO). Also the location of the most prominent kink, which comes from the heavier 
neutrino pair emission thresholds (W33 in the NO case and U12 in the IO case), can independently be 
used to extract the smallest neutrino mass value, and thus to check consistency of two experimental 
methods. 

If the spectrum is of the NO type, the measurement of the position of the kink will determine 
the value of W33 and therefore of m 3 . For the IO spectrum, the threshold U12 is very close to the 



11 



Majorana vs Dirac 
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Figure 4: Spectra from Yb 3 Pq — > 1 Sq for Dirac and Majorana neutrinos of mass mo = 20 meV. 
Dirac (black) and Majorana (red); NH (solid) and IH (dashed). 



thresholds a> 22 and Wn. The rates of emission of the pairs (z/ 2 + z/ 2 ) and iv\ + ui), however, are smaller 
approximately by the factors 10.0 and 12.7, respectively, than the rate of emission of iv\ + 1/2). Thus, 
the kink due to the iv\ + i/ 2 ) emission will be the easiest to observe. The position of the kink will 
allow to determine (mi + m 2 ) 2 and thus the absolute neutrino mass scale. If the kink due to the 
emission of {y 2 + V2) or (z/i + Ui) will also be observed, it can be used for the indivisual m 1 ,m2 
determination as well. 

The Neutrino Mass Spectrum (or Hierarchy). Once the absolute neutrino mass scale is 
determined, the distinction between the NH (NO) and IH (IO) spectra can be made by measuring the 
ratio of rates below and above the thresholds W33 and uu (or wu), respectively. We note that both 
of these measurements can be done without knowing the absolute counting rates. For m ^ 20 meV 
and NH (IH) spectrum, the ratio of the rates at u just above the W33 (wii) threshold and sufficiently 
far below the indicated thresholds, R, is given by: 

~ ; n I i I ^ — I Q33 I ^ 

NH : R(u 33 ; NH) = — ~ 2 — - = 0.70 , (38) 

2^i,j \ a ij\ 

IH : R(u u ; IH) - M 2 + 2 (h 3 | 2 + (|a 23 | 2 ) ^ Q ^ (3Q) 

\ a ij\ 

In obtaining the result f)39p in the IH case we have assumed that cj 22 and U12 are not resolved, but 
the kink due to the cun threshold could be observed. 

The latter does not corresponds to the features shown in Fig. |3] (and in the subsequent figures of 
the paper), where the kink due to the wu threshold is too small to be seen and only the kink due to 
the U\2 threshold is prominent. 

The Nature of Massive Neutrinos. The Majorana vs Dirac neutrino distinction is much more 
challenging experimentally, if not impossible, with the Yb atom. This is illustrated in Fig. HI where 
the Dirac and Majorana spectra are almost degenerate for both the NH and IH cases. The figure is 
obtained for m = 20 meV and the CPV phases set to zero, (a, (3 — 5) = (0, 0), but the conclusion is 
valid for other choices of the values of the phases as well. 
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Figure 5: The ratio -R(r) = r 72l ,(u;)/T 72! ,(u;; m; = 0) = I{uj) / I{u;mi = 0) as a function of u in 
the case of emission of Dirac and Majorana massive neutrinos having NO (left panels) or IO (right 
panels) mass spectrum with min(mj) = mo = 50; 100 meV, for e eg = 2.14 eV and four values of the 
CPV phases (a, (3 — 5) in the Majorana case. 



The difference between the emission of pairs of Dirac and Majorana neutrinos can be noticeable 
in the case of QD spectrum with mo ~ 100 meV and for values of the phases a = 0, as is illustrated in 
Fig. \5[ where we show the ratio R(T) = T^2 V {oS)/Y^2 V {u\ rrii — 0) = I{uS)/ rrii = 0) as a function 
of u. As Fig. [5] indicates, the relative difference between the Dirac and Majorana spectra can reach 
approximately 6% at values of u sufficiently far below the threshold energies Uij. For mo = 50 meV, 
this difference cannot exceed 2% (Fig. [5]). 

A lower atomic energy scale e eg > 100 meV, which is closer in value to the largest neutrino mass, 
would provide more favorable conditions for determination of the nature of massive neutrinos and 
possibly for getting information about at least some (if not all) of the CPV phases. In view of 
this we now consider a hypothetical atom X scaled down in energy by 1/5 from the real Yb, thus 
e eg ~ 0.4eV. There may or may not be good candidate atoms/molecules experimentally accessible, 
having level energy difference of order of the indicated value. Figure [6] shows comparison between 
spectra from X 3 P So for Majorana and Dirac neutrinos with m = 2 meV, for both the NH 
and IH cases. As seen in Fig. O the Majorana vs Dirac difference is bigger than 5% (10%) above 
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Majorana vs Dirac 
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Figure 6: Majorana vs Dirac comparison from X 3 P 1 So of energy difference e eg = e eg (Yb)/5 for 
rrio = 2 meV; NH in solid and IH in dashed, and the Majorana case in red, the Dirac case in black. 



the heaviest pair threshold in the NH (IH) case. The difference becomes bigger for larger values 
of the smallest neutrino mass m , making the measurement easier. This is illustrated in Fig. 
where we show again the ratio R(T) = I(uj)/I(u; = 0) as a function of u in the case of Dirac and 
Majorana pair neutrino emission for mo = 50; 100 meV and NO and IO spectra. In the Majorana 
neutrino case, the ratio R(T) is plotted for the four combinations of CP conserving values of the 
phases (a,f3 — 5) = (0,0); (0, 7r/2); (7r/2, 0); (7r/2,7r/2). There is a significant difference between the 
Majorana neutrino emission rates corresponding to (a, — 5) = (0, 0) and (vr/2, n/2). The difference 
between the emission rates of Dirac and Majorana neutrinos is largest for (a, (3 — 5) = (0,0). For 
mo = 50 (100) meV and (ot,/3 — 5) = (0,0). for instance, the rate of emission of Dirac neutrinos at 
co sufficiently smaller than CU33 in the NO case and U22 in the IO one, can be larger than the rate of 
Majorana neutrino emission by ~ 20% (70%). The Dirac and Majorana neutrino emission spectral 
rates never coincide. 

In Figs. [8] and [9] we show the spectral rate dependence on the CPV phases a and /3 — S for m = 2 
meV. Generally speaking, the CPV phase measurement is challenging, requiring a high statistics data 
acquisition. A possible exception is the case of a and IH spectrum, as shown in Fig. M, where the 
difference between the spectral rates for a = and a = n/2 can reach 10%. For the NH spectrum, 
the analogous difference is at most a few percent; observing this case requires large statistics in actual 
measurements. 

It follows from these results that one of the most critical atomic physics parameters for the 
potential of an RENP experiment to provide information on the largest number of fundamental 
neutrino physics observables of interest is the value of the energy difference e eg . Values e eg <0.4 
eV are favorable for determining the nature of massive neutrinos, and, if neutrinos are Majorana 
particles, for getting information about at least some of the leptonic CPV phases, which are the 
most difficult neutrino related observables to probe experimentally. 
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Figure 7: The same as in Fig. |5]but for e eg = 0.43 eV. 



4 Summary and conclusion 

In the present work we investigated the sensitivity of undetermined neutrino parameters (the absolute 
mass scale, the mass hierarchy pattern, the Majorana vs Dirac distinction, and the CP violating 
phases) to observables in macro-coherent RENP experiments. The easiest measurable, once RENP 
events are unambiguously identified experimentally, is the largest neutrino mass, the next easier being 
the mass hierarchy pattern (NH or IH). The Majorana vs Dirac distinction and measurement of CP 
phases are more challenging, requiring high statistics data from atoms or molecules of lower energy 
difference of order 0.5 eV. Even if measurements of these important parameters of neutrino physics 
might be demanding, a single RENP experiment may provide a systematic strategy to determine 
almost all of these parameters, thereby contributing to explore physics far beyond the standard 
theory. 

The present work points to the best atom/molecule candidate of the energy difference of less 
than O(0.5 eV) for the indicator e eg . Besides the desirable richness of detectable observables, good 
candidates for realistic RENP experiments have to be searched also from the point of easiness of 
target preparation. Investigations along these lines are in progress by a group including some of us. 
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Figure 8: CP phase dependence in NH case for the transition corresponding to Fig. El CP (a, — 5) 
assumed are (0,0) in red, (ir/2,0) in blue, and (0,7r/2) in blue dashed. 
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Figure 9: CP phase dependence in IH case for the transition corresponding to Fig. El CP (a, 13 — 5) 
assumed are (0,0) in red, (tt/2, 0) in black, and (0,7r/2) in black dashed. 
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